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Abstract 

This  paper  presents  an  advanced,  three-dimensional  model  of  the  Saturnian 
magnetosphere.  The  magnetospheric  field  of  Saturn  is  the  sum  of  three  contributions. 
The  first  two  contributions,  known  as  interior  field  sources,  are  the  inherent  planetary 
field  and  the  field  due  to  a  co-rotating  plasma  ring  around  the  planet.  The  third  field 
contribution  results  from  the  interaction  of  the  solar  wind  with  the  magnetopause  (the 
magnetospheric  boundary).  This  interaction  induces  currents  on  the  magnetopause 
surface,  which  produces  the  additional  field  contribution.  The  Cassini  Mission  is  a  joint 
venture  between  the  National  Aeronautics  and  Space  Administration  and  the  European 
Space  Agency,  designed  for  study  of  the  Saturn  system.  The  mission  planners  require  a 
magnetospheric  model  that  will  be  valid  at  the  time  of  the  Cassini  encounter  (~  June 
2004).  Previous  models  of  the  Saturnian  magnetospheric  field  either  neglected  the 
contribution  from  the  solar  wind  or  assumed  a  solar  wind  that  approached  the  planet 
parallel  to  the  plane  of  its  magnetic  equator.  These  models  do  not  depict  the  solar  wind 
orientation  that  Cassini  will  encounter  at  the  beginning  of  its  Tour,  when  the  solar  wind 
will  approach  from  the  maximum  southern  magnetic  latitude.  The  model  presented  in  this 
paper  has  been  incorporated  into  Cassini  Mission  software  to  determine  optimal  orbiter 
configuration  for  the  four- year  tour. 


Keywords:  Saturn,  Planetary  Magnetic  Fields,  Magnetosphere,  Magnetopause,  Cassini 
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I.  Introduction 

The  Cassini  Mission  is  a  joint  endeavor  between  the  National  Aeronautics  and 
Space  Administration  (NASA),  the  European  Space  Agency  (ESA),  and  the  Italian  Space 
Agency  (ASI).  The  project  is  named  in  honor  of  Jean  Dominique  Cassini  (1625-1712), 
an  Italian-born,  French  astronomer  who  discovered  the  large  division  in  Saturn’s  rings 
and  four  of  the  planet's  moons.  Scheduled  for  launch  in  October  1997  on  board  a  Titan 
IV/Centaur  rocket,  the  Cassini  orbiter  will  journey  to  Saturn  by  way  of  gravity  assists 
from  Venus,  Earth,  and  Jupiter.  The  spacecraft  will  be  injected  into  Saturn  orbit  and 
begin  its  period  of  in  situ  observations  in  June  2004,  whereupon  the  four-year  (~  60 
orbit)  data  acquisition  phase  will  begin  [Cassini  Mission  Plan,  1995]. 

The  goal  of  the  Cassini  Mission  is  to  enhance  understanding  of  the  Saturnian 
system.  To  this  end,  the  orbiter  will  focus  on  four  areas  of  scientific  investigation:  the 
planet’s  atmosphere  and  rings,  the  planetary  magnetosphere,  the  moon  Titan,  and  the 
planet’s  icy  satellites.  Cassini  will  employ  the  most  advanced  sensor  suite  ever 
assembled,  which  will  be  used  to  create  a  comprehensive  planetary  database.  The 
orbiter’ s  instrumentation  is  crucial  to  mission  success,  and  its  effectiveness  during  its  in 
situ  period  is  highly  dependent  on  the  trajectory  and  orientation  of  the  spacecraft. 

The  study  of  the  Saturnian  magnetosphere  is  the  responsibility  of  the  Cassini 
Mission’s  Magnetosphere  and  Plasma  Science  working  group  (MAPS).  To  insure 
optimal  orbiter  orientation  for  magnetospheric  studies,  MAPS  requires  a  detailed,  three- 
dimensional  model  of  the  planetary  magnetospheric  field.  The  magnetosphere  is  the 
region  around  the  planet  where  its  magnetic  field  is  the  dominant  field.  A  well-defined 
boundary  exists  between  the  planetary  magnetic  field  and  interplanetary  space.  This 
boundary  is  referred  to  as  the  magnetopause. 

Modeling  the  magnetosphere  consists  of  predicting  the  location  of  the 
magnetopause  and  describing  the  magnetic  field  it  encapsulates.  Data  collected  by 
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Pioneer  11  (PI  1),  Voyager  1  (VI),  and  Voyager  2  (V2)  revealed  that  Saturn's 
magnetosphere  is  intermediate  to  those  of  Jupiter  and  Earth.  Details  about  the  global 
magnetosphere  are  limited,  however,  because  the  trajectories  of  the  Pioneer  and  Voyager 
missions  were  limited  to  relatively  constant  magnetic  latitudes  and  only  one  (V2)  or  two 
(PI  1,  VI)  crossings  of  the  magnetic  equator  [Connerney  et  al.,  1981].  The  detailed 
magnetospheric  database  that  scientists  need  to  construct  an  advanced,  global 
magnetospheric  model  will  be  acquired  by  Cassini,  which  will  collect  data  over  a  broad 
range  of  magnetic  latitudes  and  longitudes. 

MAPS  incorporates  magnetospheric  models  into  the  Jet  Propulsion  Laboratory’s 
MAGPAC  software  package.  This  software  allows  the  mission  planning  team  to  perform 
simulated  flybys  using  different  orbital  trajectories  and  orientations.  The  results  of  these 
simulations  will  be  used  by  the  Tour  planning  team  to  determine  how  best  to  situate  the 
orbiter  for  collection  of  magnetospheric  data.  After  magnetospheric  data  collected  by  the 
orbiter  has  been  transmitted  to  Earth  and  analyzed,  the  models  used  for  Tour  planning 
will  serve  as  a  point  of  origin  for  developing  new  models  of  the  Saturnian 
magnetospheric  field. 

II.  Background 

The  Saturnian  magnetospheric  field  is  the  sum  of  three  contributions.  The  first 
two  contributions  are  referred  to  as  the  "interior  sources"  because  they  originate  within 
the  magnetosphere:  they  are  the  intrinsic  planetary  field  and  the  field  contribution  from  a 
toroidal  plasma  formation  surrounding  the  planet.  These  sources,  which  are  both  axially 
and  equatorially  symmetric,  have  been  deduced  from  Pioneer  and  Voyager  observations. 
The  third  contribution  to  the  planetary  magnetospheric  field  is  the  result  of  the  interaction 
of  the  solar  wind  with  the  magnetospheric  field  just  inside  the  magnetopause.  This 
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interaction  induces  currents  on  the  magnetopause  surface,  thus  producing  a  third  magnetic 
field  contribution. 

Analysis  of  Pioneer  and  Voyager  magnetometer  data  reveals  that  the  intrinsic 
magnetic  field  of  Saturn  is  well-represented  by  a  planetary  dipole  field  of  moment  ~  0.21 
G  •  Rs3  (1  Rs  =  60,330  km).  As  is  characteristic  of  a  dipole  field,  it  is  nearly  aligned 
(less  than  1°  offset)  with  the  planet’s  rotation  axis  [Acuna  and  Ness,  1980;  Ness  et  al., 
1981,1982].  This  alignment  is  unique  to  Saturn  when  compared  to  the  Earth  and  Jupiter, 
and  leads  to  the  near-coincidence  of  the  Saturnian  rotational  and  magnetic  planetary 
coordinate  systems.  The  dipole  field  of  Saturn  is  displayed  in  Figure  1  [Maurice  and 
Engle,  1995], 

Higher-order  models  of  the  intrinsic  field  have  been  developed  [Connemey  et  al., 
1982;  Davis  and  Smith,  1990],  but  they  diverge  from  the  dipole  model  only  at  distances 
close  to  the  planet’s  surface.  When  treating  the  entire  magnetospheric  field,  it  is  sufficient 
to  represent  the  intrinsic  planetary  field  as  a  dipole,  Bdipoie-  This  field  is  likely  the  result 
of  a  dynamo  effect,  in  which  the  metallic  hydrogen  zone  (a  region  extending  from  5000 
km  to  8000  km  from  the  center  of  the  planet)  of  the  rapidly  rotating  planet  produces  a 
magnetic  field  contribution,  dB,  given  by  the  Biot-Savart  Law  [Zeilik,  1991], 

a, 

4  7t  r2 

I  is  the  magnitude  of  the  current,  ds  is  a  differential  element  indicating  the  direction  of 
current  flow,  r  is  the  unit  vector  in  the  direction  from  the  source  point  to  the  field  point, 
and  r  is  the  distance  between  the  source  point  and  the  field  point 
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Figure  1.  The  dipole  field  of  Saturn  [Maurice  and  Engle,  1995).  The  field 
lines  lie  in  the  noon/midnight  (XZ)  plane.  Intensity  of  the  field  is  a  function 
of  the  density  of  the  field  lines.  The  X  axis  is  the  planetary  magnetic  equator 
and  in  the  positive  direction  points  to  the  Sun.  The  Z  axis  is  the  planetary 
magnetic  axis,  which  is  nearly  coincident  with  the  planetary  rotation  axis. 
The  axes  are  measured  in  units  of  Saturnian  radii  ( 1  Rs  =  60.330  km). 


After  further  analysis  of  VI  data,  Connerney  et  al.  noted  a  departure  from  the 
dipole  model  beyond  ~  8  Rs  on  the  day-side  and  outside  ~  5  Rs  on  the  night-side  of  the 
planet.  To  account  for  this  deviation  they  proposed  an  azimuthal  current  system  whose 
magnetic  field  contribution,  given  by  equation  (1),  was  supported  by  experimental  data 
[1981],  The  VI  plasma  experiment  confirmed  that  the  current  system  could  be 
represented  by  a  toroidal-shaped  plasma  region,  whose  cross  section  is  a  rectangle  with 
inner  and  outer  edges  at  8  R$  and  15.5  R$,  respectively,  and  whose  height  is  -  6  Rs 
[Connerney  et  al.,  1983]. 
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Within  the  plasma  “ring,’7  the  current  density  decreases  radially  as  1/r  from  a  value 
of  ~  0.3  x  106  A/Rs2  at  the  inner  edge  to  ~  0.15  x  106  A/Rs2  at  the  outer  edge  of  the 
formation.  This  ring,  whose  magnetic  field  contribution  will  henceforth  be  referred  to  as 
Bnng,  is  similar  in  nature  to  the  Jovian  plasma  disk  [Bode,  1994];  however,  it  is  about 
one  order  of  magnitude  smaller  in  terms  of  current  density,  contributing  ~  10  nT  to  the 
near-axis  field.  It  is  believed  that  the  plasma  is  composed  of  light  and  heavy  ions, 
probably  hydrogen,  oxygen,  or  nitrogen  [ Bridge  et  al.,  1981],  By  superposition 
principles,  Br;ng  and  Bdipoie  add  to  yield  the  total  magnetic  field  due  to  the  internal 
magnetospheric  sources.  This  field  is  displayed  in  Figure  2  [Maurice  and  Engle,  1995]. 


x(Rs) 


Figure  2.  A  view  in  the  XZ  plane  of  the  sum  of  the  dipole  field  and  the  field 
due  to  the  plasma  ring  current  (Bdipoie  +  Bnng)  [Maurice  and  Engle ,  1995]. 
The  regions  bounded  by  dotted  lines  indicate  where  the  plasma  is  contained. 
The  field  lines  originate  at  the  same  latitudes  as  those  in  Figure  1. 
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The  third  contribution  to  the  Saturnian  magnetospheric  field  arises  from  the 
interaction  of  the  solar  wind  with  the  magnetospheric  field  at  the  magnetopause.  Solar 
wind  is  a  supersonic  particle  flow  emanating  from  the  Sun  that  fills  the  interplanetary 
space  of  the  solar  system.  It  is  composed  of  approximately  equal  numbers  of  electrons 
and  protons,  such  that  there  is  no  net  current  produced  by  the  wind.  When  this  particle 
flow  reaches  the  Saturnian  magnetosphere,  the  charged  particles  of  the  flow  are  subjected 
to  a  magnetic  force,  F,  given  by 


F  =  qv  x  B.  (2) 

q  is  the  charge  on  the  particle,  v  is  the  velocity  of  the  incoming  particle,  and  B  is  the 
magnetic  field  that  the  particle  encounters. 

At  moderate  magnetic  latitudes,  this  force  redirects  protons  in  a  westward 
direction,  while  electrons  are  sent  in  an  eastward  path.  The  result  is  that  a  net  current  is 
formed  on  the  magnetopause  surface  in  the  direction  of  proton  deflection.  The  net  current 
at  each  point  on  the  surface  of  the  magnetopause  produces  an  additional  magnetospheric 
field  contribution,  BMp,  given  by  equation  (1).  This  contribution  vectorially  adds  with 
the  interior  field  sources  to  arrive  at  the  total  magnetospheric  field: 

B  =  B^ipoie  +  Bnng  +  Bmp-  (3) 


III.  Method 

The  interior  sources  of  the  planetary  magnetosphere,  Bdipoie  +  Bring*  have  an 
associated  magnetic  pressure.  Likewise,  the  change  in  momentum  of  the  solar  wind 
particles  as  they  enter  the  magnetosphere  at  the  magnetopause  exerts  a  Newtonian 
pressure.  The  balance  between  the  magnetic  pressure  of  the  interior  sources  and  the 
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incoming  solar  wind  pressure  permits  the  calculation  of  the  three-dimensional 
magnetopause  surface.  The  applicability  of  this  condition  is  discussed  in  detail  by  Beard 
[  1960]  and  Mead  and  Beard  [1964]. 

The  pressure  balance  equation  at  the  magnetopause  boundary  is 

2psVs2  cos2  X  =  •  (4) 

X  is  the  angle  between  the  incident  solar  wind  and  the  normal  to  the  magnetopause 
surface;  BT  is  the  component  of  the  net  magnetic  field  at  the  boundary,  and  is 
approximately  tangential  to  the  magnetopause;  and  ps  and  Vs  are  the  mass  density  and  the 
speed  of  the  particle  flow  [Maurice  and  Engle,  1995],  The  momentum  conservation 
required  by  equation  (4),  in  addition  to  accounting  for  all  of  the  incoming  solar  wind 
particles,  leads  to  the  following  working  equation  applicable  at  each  point  on  the 
magnetopause  surface: 


|ns  x  B  |  ±  ns  v  =  0.  (5) 

The  unit  vector  ns  is  the  normal  to  the  magnetopause  surface  given  in  spherical  polar 
coordinates  by: 


~  1  dr  n 

ns  =  Kn  r  -  4-  —  0  - 

r  ae 


3r 


r  sin  0  a<> 


(6) 


where  Kn  is  defined  as: 


K 


n 


\  dQ)  \r  sin  0  dty) 


(7) 
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The  unit  vector  v  delineates  the  incoming  solar  wind  direction.  B  is  the  total  magnetic 
field  vector,  normalized  to  unity  at  the  subsolar  point,  which  is  where  the  incoming  solar 
wind  is  perpendicular  to  the  magnetopause  surface.  In  the  first  calculation  of  the 
magnetopause  surface,  B  consists  of  only  Bdipoie  +  B^g  because  the  contribution  from 
the  surface  currents  depends  on  that  surface’s  curvature,  which  is  not  yet  calculated 
[Mead  and  Beard,  1964]. 

Following  construction  of  the  trial  surface,  the  first  set  of  normalized 
magnetopause  surface  current  elements  can  be  computed  according  to: 

J=£sxB.  (8) 


Substituting  equations  (6)  and  (7)  into  equation  (8),  one  arrives  at  an  alternate  expression 
for  J: 

\r  sin  0  3(f)  30  /  \  r  sin  0  3<J)  /  \  30  / 


After  these  current  elements  have  been  computed  for  the  entire  magnetopause  surface, 
they  may  be  used  to  calculate  a  corresponding  magnetic  field  contribution,  Bsc-  This 
contribution  is  given  by  the  Biot-Savart  relationship  in  the  following  form: 

Bosc=l /Jxirr^dS'  (10) 

By  substitution  of  equation  (8)  into  equation  (10),  the  following  alternate  expression  is 


derived: 
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B„,  sc  =  f  f  (ns  x  Bn .  i)  x  dS’  (11) 

The  procedure  (5)  -  (11)  is  performed  iteratively,  incorporating  new  magnetopause 
surfaces  and  the  corresponding  current  elements  until  the  model  becomes  self-consistent. 
That  is, 


Bn.  sc  =  Bn-i.  sc 


(12) 


and,  consequently. 


Bn  =  Bn-!.  (13) 

The  currents  on  the  magnetopause  have  a  corresponding  magnetic  field 
contribution,  Bmp-  At  points  inside  the  magnetosphere,  the  magnetic  field  may  be 
represented  as  the  gradient  of  a  scalar  magnetic  potential,  given  in  spherical  polar 
coordinates  by 


=  -V<I>Mp(r,0,<t>). 


(14) 


A  possible  expansion  of  <J>Mp(rA<t>)  is 


Umax  n 

<t>MP(r,e,<|>)  =  X  141  X  °nm  P{?(cos  0)  [cos  (m<(>)]. 

n  =  1  m  =  0 


(15) 
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The  coefficients  Gnm  are  determined  by  using  the  calculated  values  of  BMP  at  points 
inside  the  magnetopause.  This  allows  one  to  obtain  the  total  magnetic  field  of  the  planet, 
given  by 


B  =  B 


dipole 


+  Bring  +  B 


MP- 


(16) 


IV.  Existing  models 

Data  collected  by  the  magnetometers  of  PI  1,  VI,  and  V2  has  been  used  as  the 
basis  for  the  construction  of  several  models  of  the  Saturnian  magnetospheric  field.  Acuna 
et  al.  [1980]  present  a  field  model  based  on  analysis  of  Pll  data  [Acuna  and  Ness, 
1979],  and  Connerney  et  al.  [1981]  present  a  model  constructed  from  analysis  of  VI  data 
[Ness  et  al.,  1981].  Each  of  these  models  was  based  on  two  assumptions:  a  dipole 
representation  of  the  inherent  planetary  field  and  the  alignment  of  the  planet’s  rotation  and 
magnetic  axes. 

The  planetary  rotational  and  magnetic  axes  are  not  coincident  with  the  magnetic 
axis  offset  by  ~  0.75°  [New  etal.,  1981].  The  Zonal  harmonic  (Z3)  model,  presented  by 
Connerney  et  al.  [1982],  and  the  Saturn  Pioneer  Voyager  (SPV)  model,  constructed  by 
Davis  and  Smith  [1990],  attempt  to  account  for  the  field  contributions  from  the  slight 
variation  in  alignment  of  the  planetary  rotational  and  magnetic  axes.  The  variations 
between  these  models  and  the  dipole-like  models  are  minor,  especially  at  larger  distances 
from  the  planet.  Additionally,  these  models  do  not  account  for  the  magnetic  field 
contribution  associated  with  the  solar  wind  interaction  at  the  magnetopause. 

The  PI  1,  VI,  and  V2  encounters  with  the  Saturnian  magnetosphere  were  brief- 
on  the  order  of  hours— thereby  limiting  the  potential  to  collect  data.  In  addition  to  time 
constraints,  magnetometer  data  collected  by  the  probes  was  limited  to  small  magnetic 
latitudes  because  of  mission  trajectories.  Thus  V 2  crossed  the  magnetic  equatorial  plane 
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only  once,  while  Pll  and  VI  each  experienced  only  two  crossings.  The  short  in  situ 
periods  and  the  mission  trajectories  are  the  two  factors  that  have  limited  the  collection  of 
Saturnian  magnetospheric  data.  This,  in  turn,  limits  the  aforementioned  models  to  a  small 
range  of  magnetic  latitudes  [Connemey  et  al.,  1981]. 

V.  The  Maurice/Engle  idealized  model 

To  extend  the  scope  of  knowledge  of  the  Saturnian  magnetosphere,  Maurice  and 
Engle  constructed  the  first  advanced  three-dimensional  model  of  the  Saturnian 
magnetosphere  [1995],  This  model  represents  the  Saturnian  magnetospheric  field  based 
on  the  method  outlined  in  Section  III  (although  the  method  has  been  used  to  model  the 
fields  of  Earth  and  Jupiter)  [Beard,  1960;  Engle  and  Beard,  1980]. 

The  model  constructed  by  Maurice  and  Engle  (subsequently  referred  to  as  ME95) 
incorporated  interior  field  sources  due  to  the  intrinsic  planetary  field  and  the  ring  current. 
Existing  models  attempted  to  represent  these  fields  as  well;  however,  this  model  differs 
from  previous  efforts  because  it  took  into  account  the  contribution  from  the  magnetopause 
surface  currents.  The  result  is  that  the  physically  reliable  range  of  the  model  was 
extended  to  nearly  all  magnetic  latitudes  and  longitudes  within  the  Saturnian 
magnetosphere. 

The  Pioneer  and  Voyager  missions  encountered  a  large  variability  in  the  magnetic 
environment  around  Saturn.  To  account  for  this,  ME95  modeled  the  magnetospheric 
field  based  on  six  parameters  that  describe  the  magnetosphere:  one  for  the  planetary  field, 
four  for  the  ring  current,  and  one  for  the  strength  of  the  incoming  solar  wind.  In  their 
model,  Maurice  and  Engle  represented  the  planetary  field  by  a  dipole  with  moment  of  ~ 
0.21  G  •  Rs3.  The  axially  symmetric  ring  current  proposed  by  Connerney  et  al.  [1983] 
was  incorporated,  with  a  =  8  Rs  and  b  =  15.5  R$.  The  average  current  strength  was  set 


at  l0  =  2.9  x  106  A/Rs,  which  is  independent  of  its  distance  from  the  magnetic  equator, 
|  z  |  ^  D  =  3  Rs-  The  subsolar  point  was  chosen  to  be  R0  =  20  Rs. 

The  parameters  a,  b,  and  D  fix  the  dimensions  of  the  toroidal  plasma  formation  in 
Saturn's  inner  magnetosphere.  With  the  corresponding  current  density  I0,  the  total 
current  of  the  ring  may  be  calculated  and  its  contribution  to  the  total  magnetospheric  field 
deduced.  The  subsolar  point  R0  is  a  parameter  that  represents  the  relative  strength  of  the 
solar  wind.  Depending  on  the  pressure  of  the  wind,  the  location  of  the  subsolar  point 
will  vary.  The  model  presented  by  Maurice  and  Engle  [  1995]  is  displayed  in  Figure  3. 
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Figure  3.  The  noon/midnight  view  of  the  idealized  magnetospheric  model  constructed 

by  Maurice  and  Engle  [1995].  The  dotted  line  is  the  magnetopause  surface.  The  solar 

wind  is  coming  from  the  right,  parallel  to  the  X  axis.  This  introduces  symmetry  about 
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ME95  is  referred  to  as  an  idealized  model  because  of  the  direction  of  the  incoming 
solar  wind,  which  approaches  parallel  to  the  plane  of  the  magnetic  equator.  This 
introduces  a  natural  symmetry  between  the  northem/southem  hemispheres  and  the 
eastern/westem  hemispheres  when  calculating  the  magnetopause  boundary.  The 
calculation  of  the  magnetopause  surface  is,  therefore,  reduced  to  one  quarter  of  its  total 
surface  area.  From  that,  the  other  three  regions  could  be  obtained  immediately  by 
applying  symmetry  considerations.  While  the  solar  wind  approached  in  this  manner  at 
the  time  of  the  Pioneer  and  Voyager  encounters  with  Saturn,  this  orientation  will  not  be 
prevalent  when  Cassini  arrives  at  the  planet  in  June  2004. 

VI.  The  solar  wind  variation 

The  orbital  plane  containing  both  Satum  and  the  Sun  is  referred  to  as  the  ecliptic 
plane.  The  rotation  axis  (axis  about  which  the  planet  spins)  of  Satum  makes  an  angle  of 
26.7°  with  the  normal  to  the  ecliptic.  This  orientation  is  constant  in  space  throughout  the 
planet’s  orbit,  leading  to  a  cyclical  variation  of  the  rotation  axis  orientation  with  respect  to 
the  Sun  over  the  course  of  a  Satum  year  (29.46  Earth  years).  The  rings  of  Satum  lie  in 
the  planet’s  equatorial  plane,  to  which  the  rotation  axis  is  normal;  therefore,  the  cyclical 
variation  of  the  orientation  of  the  rings  with  respect  to  the  Sun  corresponds  to  the 
variation  of  the  rotation  axis  orientation  during  a  Satum  year. 

The  change  of  Saturn’s  rotation  axis  orientation  with  respect  to  the  Sun  is 
analogous  to  the  cause  of  seasons  on  Earth.  Consider  the  First  Saturnian  season 
beginning  when  the  rotation  axis  points  toward  the  Sun.  In  this  configuration,  the  rings 
of  the  planet  will  dip  down  when  viewed  from  the  Earth  (or  the  Sun),  making  an  angle  of 
-26.7°  with  the  ecliptic  plane.  A  quarter  cycle  later,  the  axis  points  neither  at  the  Sun  nor 
away  from  the  Sun.  This  neutral  orientation  means  that  the  rings  of  the  planet  may  be 
viewed  edge-on  from  Earth.  Next,  the  planet’s  rotation  axis  points  away  from  the  normal 
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to  planet-sun  line.  During  that  era,  the  rings  will  appear  tilted  upward  when  viewed  from 
the  Earth.  After  the  planet  has  proceeded  through  the  last  quarter  of  its  orbit,  the  axis  will 
once  again  be  perpendicular  to  the  sun-planet  line  and  the  rings  may  be  viewed  edge-on. 

Based  on  Pioneer  and  Voyager  data,  it  has  been  determined  that  the  rotational  and 
magnetic  coordinate  systems  of  Saturn  are  nearly  coincidental,  with  the  planetary 
magnetic  axis  offset  from  the  planetary  rotation  axis  by  less  than  1°.  As  a  result,  the 
magnetic  coordinate  frame  of  the  planet  completes  one  full  cycle  with  respect  to  the  Sun  in 
the  same  manner  as  does  the  rotational  frame,  namely,  once  every  Saturnian  year.  This  is 
a  trait  unique  among  planetary  magnetospheres.  The  axis  of  Jupiter,  for  example,  makes 
an  angle  of  10°  relative  to  the  planetary  rotation  axis.  Thus,  the  Jovian  magnetosphere 
precesses  about  the  planetary  rotation  axis  on  a  daily  (1  Jovian  day  =  9h  55m)  basis 
[. Bode ,  1994], 

The  configuration  in  which  the  rings  of  Saturn  can  be  viewed  edge-on  occurs 
twice  during  a  Saturnian  year.  This  was  the  orientation  during  the  Pioneer/Voyager 
epoch  (~  1980)  and  again  in  1995,  half  an  orbital  cycle  later.  The  ME95  model 
incorporates  this  orientation,  which  corresponded  to  data  taken  during  the 
Pioneer/Voyager  time  frame.  However,  when  Cassini  arrives  at  the  planet  (~  June  2004), 
the  planetary  orientation  will  be  such  that  the  rotation  axis  is  pointed  away  from  the  Sun 
(the  rings  are  tilted  up  when  viewed  from  the  Earth).  This  presents  a  tour  planning 
dilemma  for  the  Cassini  Mission. 

The  team  planning  Cassini’s  Tour  requires  a  model  of  the  magnetosphere  of 
Saturn  that  can  be  anticipated  at  the  beginning  of  the  probe’s  in  situ  period.  The  changes 
in  the  Saturnian  magnetosphere  are  related  to  the  variation  of  the  planet’s  position  in  orbit 
around  the  Sun.  In  order  to  anticipate  this  condition,  one  must  anticipate  the  angle  of 
incidence  that  the  solar  wind  will  have  upon  the  Saturnian  magnetopause,  since  it  is  a 
significant  factor  affecting  the  shape  of  the  planetary  magnetosphere. 
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While  the  solar  wind  does  not  emanate  radially  from  the  Sun,  the  distance  that  it 
travels  before  reaching  Saturn  insures  that  vectors  representing  the  solar  wind  velocity 
will  be  parallel.  Based  on  this  fact,  it  is  convenient  to  introduce  the  parameter  X.  When 
the  rings  of  Saturn  are  viewed  edge-on  from  the  Earth,  the  solar  wind  approaches  parallel 
to  the  plane  of  the  planetary  magnetic  equator.  This  corresponds  to  an  angle  of  X  =  0°. 
When  the  Cassini  probe  arrives  at  the  planet,  the  rings  will  be  near  their  maximum 
upward  tilt  as  viewed  from  the  Earth.  Thus,  the  angle  X  will  take  the  same  value  as  the 
angle  between  the  plane  of  the  rings  and  the  ecliptic  plane  .  As  discussed  previously,  the 
solar  wind  will  be  incident  at  the  largest  possible  value  south  of  the  magnetic  equator, 
namely,  at  X  =  -  26.7°. 

VII.  Constructing  the  X  =  -  26.7°  model 

This  paper  presents  a  magnetospheric  model  utilizing  the  condition  X  =  -  26.7°. 
The  model,  which  anticipates  the  conditions  at  Saturn  just  prior  to  the  arrival  of  Cassini, 
may  be  used  by  MAPS  scientists  at  the  Jet  Propulsion  Laboratory  (JPL)  in  Cassini 
Mission  simulation  software,  called  MAGPAC.  After  MAPS  converts  the  model  into 
computer  code,  it  is  loaded  into  MAGPAC  to  perform  simulated  flybys  of  the  planet. 
MAPS  will  determine  the  optimal  orbiter  configuration  for  magnetospheric  data 
acquisition  based  on  the  MAGPAC  simulations.  Once  this  configuration  is  chosen,  it  will 
be  considered  with  the  requests  of  other  mission  scientists  to  fix  the  exact  Tour  of  the 
spacecraft. 

Construction  of  the  X  =  -  26.7°  model  presented  here  began  in  a  similar  fashion  to 
ME95.  The  magnetospheric  environment  of  Saturn  was  again  characterized  by  the  six 
minimum  parameters.  The  planetary  field  was  represented  as  a  dipole  moment  -  0.21 
G  •  Rs3.  The  plasma  current  formation  extended  from  a  =  8  Rs  to  b  =  15.5  Rs,  with  a 
current  density  parameter  of  2.9  x  106  A/Rs.  The  height  of  the  plasma  ring  above  the 
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equatorial  plane  was  |  z  |  <  D  =  2.8  Rs.  The  subsolar  point  was  chosen  to  be  R„  =  24 

Rs 

These  parameters  were  fixed  based  on  the  parameters  used  in  the  ME95  model. 
Slight  modifications  were  made,  however,  to  scale  the  parameters  to  a  different  subsolar 
point.  This  was  necessary  because  the  subsolar  point,  while  being  indicative  of  the 
strength  of  the  solar  wind,  is  also  representative  of  the  value  of  the  cutrent  circulating  in 
the  plasma  ring.  Based  on  the  best  fit  to  VI  data,  the  subsolar  distance  *„  =  24  Rs  was 
chosen.  This  subsolar  distance  is  consistent  with  the  observed  variauons  of  the  subsolar 
point,  according  to  VI  observations. 

Construction  of  the  model  first  focused  on  building  a  trial  magnetopause  surface. 
Using  a  5°  x  5°  grid  as  a  framework  for  numerical  calculations,  construction  began  using 
FORTRAN  code  on  a  Macintosh  Quadra  800.  First  the  noon  meridian,  which  lies  in  the 
XZ  plane  in  the  planetary  magnetic  frame  of  reference,  was  calculated.  Following  this, 
the  midnight  meridian  was  computed.  These  meridia,  because  of  the  near-coincidence  of 
the  magnetic  and  rotation  axes,  inherently  yielded  a  plausible  surface  meridia  framework 
from  which  computation  of  the  remainder  of  the  surface  was  accomplished. 

Construction  of  the  ME95  model  was  required  over  only  one  quadrant  of  the 
magnetopause  surface.  The  surface  could  then  be  extended  to  the  other  three  quadrants 
by  due  to  its  four-fold  symmetry.  However,  for  the  X.  =  -  26.7°  surface,  the  only 
symmetry  present  is  that  between  the  eastern  and  western  hemispheres.  As  a  result, 
calculations  were  performed  over  a  surface  area  twice  as  large  as  the  ME95  model.  The 
breakdown  of  symmetry  between  the  northern  and  southern  hemispheres  introduced 
programming  challenges  not  encountered  in  the  idealized  situation. 

While  constructing  the  trial  magnetopause  surface,  the  computer  calculated  the 
difference  between  the |  ns  x  B  |  and  the  |ns  •  v|  terms  of  equation  (5).  The  solution, 
which  is  a  distance  r  from  the  planet,  was  chosen  based  on  where  this  difference  was 
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minimized.  In  attempting  to  find  the  surface  point,  the  program  first  looked  for  a  solution 
based  on  points  in  the  neighboring  region.  If  this  failed,  or  no  points  had  been 
calculated,  the  program  set  a  maximum  and  minimum  value  for  r .  The  lower  limit  was 
based  on  the  calculation  of  the  previous  point,  while  the  upper  limit  involved 
extrapolating  a  line  perpendicular  to  the  normal  at  the  magnetopause  surface  of  the 
previously  calculated  point. 

The  computer  could  be  left  to  calculate  the  magnetopause  surface  on  its  own; 
however,  the  vast  majority  of  its  results  would  not  be  physically  plausible.  Scientists 
have  attempted  to  automate  the  process,  but  due  to  the  number  of  iterations  required  and 
the  multitude  of  possible  solutions  involved  in  calculating  a  self-consistent  magnetopause 
surface,  however,  these  attempts  have  failed.  The  process  requires  the  physical  intuition 
of  the  physicist  in  order  to  choose  numerical  solutions  that  coincide  with  the  current 
knowledge  of  planetary  magnetospheres  and  the  relatively  vague  boundary  conditions, 
which  introduce  families  of  numerically  valid  but  physically  unacceptable  solutions. 

When  the  model  was  not  near  a  state  of  self-consistency,  solutions  that  yielded  a 
more  physically  feasible  surface  were  not  selected  by  the  program,  which  chooses  the 
best  numerical  solution.  To  compensate  for  this,  it  was  necessary  to  individually 
investigate  solutions  that  posed  physically  unfeasible  situations.  This  was  done  by 
inspecting  the  program  outputs  in  an  attempt  to  find  local  minima  in  the  difference 
between  the  two  terms  of  equation  (5)  that  the  computer  may  have  overlooked.  This  was 
a  crucial  step  in  computing  the  first  several  iterations  of  the  magnetopause  surface. 

VIII.  Results  of  the  X  =  -  26.7°  model 

Self-consistency  of  the  X  =  -  26.7°  model  was  achieved  after  ten  iterations.  The 
magnetopause  surface  that  was  calculated  is  displayed  in  Figure  4.  The  asymmetry 
effects  of  the  solar  wind  are  apparent  when  the  X  =  -  26.7°  magnetopause  surface  is 
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compared  to  the  ideal  surface  of  ME95.  The  incoming  solar  wind  pressure  tends  to 
compress  the  region  at  the  front  of  the  magnetopause  just  below  the  day-side  magnetic 
equatorial  plane.  Additionally,  the  solar  wind  pushes  the  magnetotail  far  above  the 
magnetic  equator  of  the  planet  on  the  dusk-side. 


After  computing  a  self-consistent  magnetopause  surface,  an  expansion  in 
associated  Legendre  functions  yielded  the  coefficients  Gnm  that  reproduce  the  planetary 
field  lines  according  to  equation  (15),  outlined  in  Section  III.  Appendix  B  lists  the 
coefficients  Gm  necessary  to  reproduce  the  magnetospheric  field.  This  field  is  displayed 
in  Figure  5. 
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Figure  5.  Field  lines  of  the  Saturnian  magnetospheric  field  in  the 
noon/midnight  plane  when  X.  =  -  26.7°.  The  solar  wind  approaches  from  the 
right  and  makes  an  angle  of  -  26.7°  with  the  X  axis. 


When  compared  to  the  ME95  idealized  model,  the  effect  of  the  solar  wind 
asymmetry  upon  the  shape  of  the  magnetosphere  is  evident,  particularly  in  the  "cusp" 
regions  of  the  field.  The  "cusp"  is  the  region  of  the  magnetosphere  where  the  field  lines 
cease  to  reconnect  to  the  planet  and,  instead,  sweep  back  to  form  the  magnetotail.  While 
the  interior  field  near  the  magnetic  equator  is  similar  to  that  of  the  ME95  model,  the  cusp 
regions  are  quite  different  in  terms  of  relative  north  and  south  field  intensities  and 
latitudes  at  which  they  occur.  As  the  cusp  regions  are  generally  associated  with  planetary 
auroral  phenomena,  where  charged  particles  spiral  along  magnetic  field  lines  toward  the 
planet  surface,  asymmetrical  auroral  phenomena  should  be  anticipated  in  future 
observations. 
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Consideration  of  the  field  lines  close  to  the  planet’s  surface  shows  that  the  solar 
wind  apparently  both  compresses  the  day-side  field  and  pushes  the  night-side  field  lines 
upward.  This  compression  may  be  seen  in  Figure  6.  These  field  lines  originate  at  the 
same  magnetic  latitude  (S  80°)  in  the  region  of  the  southern  cusp  and  are  calculated 
around  half  of  the  planetary  surface  in  30°  increments.  The  day-side  field  lines  are 
clearly  compressed  by  the  solar  wind,  while  the  night-side  lines  appear  relatively 
extended  and  pushed  upward.  The  effect  is  most  noticeable  on  the  fourth  field  line  from 
the  right,  which  originates  in  the  YZ  plane.  As  the  line  travels  further  from  the  planet's 
surface  it  is  pushed  out  of  the  YZ  plane  by  the  solar  wind  pressure. 


Figure  6.  Field  lines  originating  at  a  magnetic  latitude  of  S  80°.  The  lines  begin  in  the  noon  (positive 
X)  plane,  continuing  in  30°  increments  to  the  midnight  (negative  X)  plane. 
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IX.  Validating  the  X  =  -  26.7°  model 

The  in  situ  data  that  has  been  collected  on  the  Saturnian  magnetospheric  field  to 
date  corresponds  to  measurements  that  were  taken  during  the  epoch  of  Pioneer  and 
Voyager.  Thus,  because  the  solar  wind  was  configured  such  that  X  ~  0°,  no  direct 
observations  are  available  to  confirm  the  predictions  of  the  X  =  -  26.7°  model.  There  are, 
however,  two  ways  that  one  may  infer  the  accuracy  of  the  model  presented.  One 
involves  the  successful  application  of  the  method  used  here  to  calculate  the 
magnetospheric  fields  of  Earth  and  Jupiter,  while  the  other  involves  comparison  of  the 
magnetopause  surfaces  of  Saturn  as  X  varies  from  -  26.7°  to  0°. 

The  method  outlined  in  Section  HI  has  been  applied  to  both  Earth  [Beard,  1960; 
Mead  and  Beard,  1964]  and  Jupiter  [Engle  and  Beard,  1980;  Engle,  1991],  In  each  case 
the  field  model  was  consistent  with  magnetospheric  field  data  that  has  been  collected. 
The  Jovian  model  constructed  using  this  method  helped  successfully  explain  the  behavior 
of  ionized  effluvia  of  comet  fragments  from  Shoemaker-Levy  9  prior  to  their  impact  with 
Jupiter  in  July  1994  [Prange  et  al.,  1995]. 

Comparison  of  the  ME95  model  with  VI  observations  also  lends  validity  to  the 
model  presented  here.  Through  addition  of  the  magnetopause  surface  current 
contribution  to  the  planetary  magnetic  field,  Maurice  and  Engle  were  able  to  greatly 
improve  the  fit  between  the  observed  data  and  the  predicted  values  for  the  behavior  of  the 
magnetosphere  in  the  outer  region  of  the  Saturnian  magnetosphere  (17  Rs  <  r  <  24  Rs) 
[1995].  This  region  will  be  the  focus  of  much  of  Cassini’s  magnetospheric  studies. 
[Cassini  Mission  Plan,  1995]. 

Given  the  applicability  of  the  ME95  model,  a  natural  progression  of 
magnetopause  surfaces  between  the  X  =  -  26.7°  and  the  X  =  0°  models  would  be  a 
reasonable  indication  of  the  feasibility  of  the  model  presented  in  this  paper.  This  was 
done,  beginning  with  the  self-consistent,  X  =  -26.7°  magnetopause  surface.  Then  the 
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value  of  X  was  then  varied  in  an  ascending  fashion  to  X  =  -  20°,  -  15°,  -  10°,  and  -  5°. 
Trial  magnetopause  surfaces  were  constructed  for  each  of  these  values  of  X.  The  surfaces 
were  not  all  calculated  to  the  point  of  self-consistency;  therefore,  they  are  merely  an 
attempt  to  observe  the  behavior  of  the  magnetopause  as  the  angle  of  incidence  of  the 
incoming  solar  wind  increases  from  its  minimum  to  the  ideal  case.  The  results  may  be 
seen  in  Figure  7. 


Figure  7.  The  magnetopause  surface  calculations  for  X  =  -  26.7°,  -  20°,  -  15°, 
-  10°,  -  5°,  and  0°.  For  the  X  =  -  26.7°  model  to  be  physically  feasible,  the 
magnetopause  surfaces  must  gradually  transition  between  the  extreme  states, 
shown  in  the  solid  lines. 


Figure  7  indicates  a  smooth  transition  between  the  extreme  and  ideal  solar  wind 
angles.  As  X  is  varied  to  the  ME95  condition,  the  magnetopause  surfaces  gradually  shift 
between  the  ideal  and  the  X  =  -  26.7°  configuration.  This  can  best  be  seen  by  the  varying 
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location  of  the  cusp  and  the  shift  of  the  magnetotail  toward  its  symmetrical  configuration 
about  the  X  axis.  Each  of  the  transition  states  presented  is  more  "idealized"  than  the 
previous,  an  indication  of  the  physical  feasibility  of  the  surfaces.  Thus,  it  appears  that  the 
X  =  -  26.7°  model  is  physically  feasible. 

X.  Future  applications 

The  model  presented  is  the  first  in  a  series  of  steps  to  produce  a  time  dependent 
model  that  will  effectively  predict  the  magnetospheric  configuration  at  any  point  in  the 
planet's  orbit  (that  is,  for  any  angle  of  incidence  of  the  incoming  solar  wind).  First,  the 
development  of  self-consistent  surfaces  from  the  trial  surfaces  displayed  in  Figure  7 
would  yield  field  models  and  their  associated  expansion  coefficients,  Gnn,.  Once  these 
expansion  coefficients  were  known,  they  could  be  represented  as  functions  of  time  that 
would  successfully  interpolate  between  values,  predicting  the  coefficients  of  the  various 
field  models.  This  would  eliminate  the  static  nature  of  the  model,  making  it  instead  a 
function  of  time  capable  of  predicting  the  magnetospheric  field  for  any  angle  X  based  on 

the  planet’s  position  in  its  orbital  cycle. 

The  tilted  solar  wind  magnetopause  surfaces  displayed  in  Figure  7  are  also  being 
used  to  parameterize  various  configurations  of  the  Saturnian  magnetopause  [Maurice  et 
al.,  1996].  This  is  done  by  fitting  ellipses  to  the  surfaces  in  the  YZ  plane.  These  fits 
allow  one  to  predict  a  value  for  r  for  any  point  on  the  magnetopause  surface.  This  is 
significant  because  the  method  of  the  model  presented  limits  knowledge  of  r  to  the  nearest 
point  on  the  5°  x  5°  grid  that  has  been  calculated.  This  functional  representation  of  the 
magnetopause  surface  greatly  extends  our  knowledge  of  the  behavior  of  the 
magnetospheric  boundary. 

Advanced  applications  in  the  study  of  auroral  phenomena,  similar  to  those  of 
Jupiter,  are  facilitated  by  the  presentation  of  the  X  =  -  26.7°  magnetospheric  model.  As 
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previously  stated,  the  interactions  responsible  for  auroral  phenomena  are  believed  to  be 
the  result  of  charged  particles  trapped  along  field  lines  in  the  inner  magnetosphere  of  the 
planet.  In  order  to  follow  the  regions  in  which  charged  particles  may  be  accelerated  by 
the  field,  it  is  imperative  to  know  where  on  the  surface  the  field  lines  originate  and 
terminate.  This  permits  the  prediction  of  the  location  of  auroral  phenomena  so  that  further 
research  may  take  place. 

XI.  Conclusion 

An  ideal,  three  dimensional  model  of  the  Saturnian  magnetosphere  has  been 
presented.  Previous  models  either  ignored  the  solar  wind  contribution  to  the 
magnetospheric  field  or  examined  the  solar  wind  configuration  that  was  prevalent  during 
the  Pioneer/Voyager  epoch  (1979-1981).  The  model  presented  differs  from  previous 
efforts  in  that  it  utilizes  an  incoming  solar  wind  angle  of  \  =  -  26.7°.  This  is  the 

anticipated  configuration  of  the  solar  wind  at  the  time  of  injection  of  Cassini  into  orbit  at 
Saturn  (-  June  2004). 

The  model  is  required  for  use  with  the  Cassini  Tour  planning  team’s  MAGPAC 
software  package,  which  allows  the  team  to  simulate  flybys  of  the  orbiter  through  its 
Tour  utilizing  various  orientations.  This  will  permit  scientists  to  determine  the  optimal 
spacecraft  configuration.  The  motivation  for  the  optimum  orientation  is  for  the  most 
efficient  and  productive  collection  of  data  that  the  orbiter  can  accomplish  during  its  four- 
year  tour.  Currently,  the  model  presented  is  being  used  by  the  Magnetosphere  and 
Plasma  Science  working  group  of  the  Cassini  Mission  in  tour  simulation. 

While  the  model  is  based  on  primarily  theoretical  calculations,  the  same  process 
has  been  used  successfully  to  model  the  magnetospheric  fields  of  both  the  Earth  and 
Jupiter.  Additionally,  trial  magnetopause  surfaces,  obtained  by  varying  the  angle  of 
incidence  of  the  incoming  solar  wind  from  the  maximum  southern  approach  to  the  ideal 
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configuration  of  the  ME95  model,  implies  that  the  model  is  physically  feasible.  This 
indicates  that  the  model  successfully  models  the  magnetospheric  field  for  the  conditions 
that  will  exist  just  before  the  Cassini  encounter. 

Continuation  of  the  work  on  the  model  presented  will  lead  to  magnetospheric 
models  whose  expansion  coefficients  can  be  replaced  by  functions  in  time  that  will 
successfully  predict  static  field  models.  This  will  permit  calculation  of  the  field  for  any 
orientation  of  the  incoming  solar  wind  as  a  function  of  time.  The  model  is  also  being 
used  in  work  attempting  to  parameterize  Saturnian  magnetopause  surfaces.  This  will 
allow  the  distance  of  the  magnetopause  to  be  calculated  at  any  point  over  its  surface, 
rather  than  the  5°  x  5°  grid  that  is  currently  possible. 
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APPENDIX  A:  FREQUENTLY  USED  SYMBOLS 

Symbols  used  in  the  text  are  listed  below,  in  order  of  appearance. 


Symbol 


Definition 


G 

Rs 


Bdipole 

dB 

Ho 

I 

ds 


r 

r 

®ring 

F 

q 

V 

Bmp 

Ps 

Vs 

X 

By 


ru 


v 


gauss;  unit  of  measure  of  magnetic  field;  1  T  =  104  G 

Saturnian  radius;  -  60,330  km 

magnetic  field  contribution  from  planetary  dipole  field 

differential  magnetic  field  element 

permeability  of  free  space;  4rt  x  10* 7  N/A2 

current  magnitude 

differential  directional  current  element 

unit  vector  in  direction  from  source  point  to  field  point 

distance  from  source  point  to  field  point 

magnetic  field  contribution  from  co-rotating  plasma  ring 

force 

charge 

velocity 

magnetic  field  of  the  magnetopause 
mass  density  of  solar  wind  particles 
speed  of  solar  wind  particles 

angle  between  solar  wind  and  normal  to  the  magnetopause 
tangential  component  of  magnetic  field  at  the  magnetopause 
unit  vector  normal  to  the  magnetopause 
unit  vector  in  direction  of  solar  wind 
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Symbol 

Definition 

J 

current  density 

Bsc 

magnetic  field  contribution  from  magnetopause  surface  currents 

Bn 

n*  iteration  of  the  magnetic  field 

r 

distance  vector  from  origin  to  source  point 

T' 

distance  vector  from  origin  to  field  point 

dS' 

differential  distance  element 

V 

gradient 

^MP 

scalar  potential  corresponding  to  Bmp 

Gnm 

expansion  coefficients  for  the  total  magnetospheric  field 

P?  (cos  0) 

Associated  Legendre  polynomials 

a 

inner  radius  of  plasma  region,  as  measured  from  planet’s  center 

b 

outer  radius  of  plasma  region,  as  measured  from  planet’s  center 

D 

height  of  plasma  region  above  equatorial  plane 

lo 

current  density  parameter 

Ro 

subsolar  point,  as  measured  from  planet’s  center 

APPENDIX  B:  THE  COEFFICIENTS  Gnm 
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The  expansion  coefficients  used  in  equation  (15)  to  produce  the  magnetospheric  field 
line  models  for  the  conditions  X  =  0°  [Maurice  and  Engle,  1995]  and  X  =  -26.7°. 


G  (n,m) 


G  (1,0) 
G  (1,1) 
G  (2,0) 
G  (2,1) 
G  (2,2) 
G  (3,0) 
G  (3,1) 
G  (3,2) 
G  (3,3) 
G  (4,0) 
G  (4,1) 
G  (4,2) 
G  (4,3) 
G  (4,4) 
G  (5,0) 

G  (5,1) 

G  (5,2) 

G  (5,3) 

G  (5,4) 

G  (5,5) 


X=0° 

6.58E-01 

0 

0 

3.34E-01 

0 

2.23E-02 

0 

6.85E-02 

0 

0 

7.19E-02 

0 

8.21E-03 

0 

-8.70E-03 

0 

3.63E-02 

0 

-2.26E-05 


X  =  -26.70 


6.61E-01 

-5.34E-02 

-2.10E-01 

2.74E-01 

-7.86E-03 

7.75E-02 

-8.99E-02 

6.27E-02 

-7.03E-03 

-4.26E-02 

6.67E-02 

-1.22E-02 

1.89E-02 

-7.28E-03 

3.65E-03 

-3.89E-02 

2.34E-02 

1.18E-02 

1.29E-02 


0 


-3.97E-03 


G  (n,m) 


X  =  0° 


G  (6,0) 
G  (6,1) 
G  (6,2) 
G  (6,3) 
G  (6,4) 
G  (6,5) 
G  (6,6) 
G  (7,0) 
G  (7,1) 
G  (7,2) 
G  (7,3) 
G  (7,4) 
G  (7,5) 
G  (7,6) 
G  (7,7) 
G  (8,0) 
G  (8,1) 
G  (8,2) 
G  (8,3) 
G  (8,4) 
G  (8,5) 
G  (8,6) 


X  =  -26.7° 


0 

1.82E-02 

0 

1.35E-02 

0 

-5.54E-05 

0 

-1.26E-02 

0 

1.59E-02 

0 

5.03E-03 

0 

-3.39E-04 

0 

0 

-5.19E-05 

0 

5.48E-03 

0 

1.99E-03 

0 


2.77E-03 
6.31E-03 
-1.89E-02 
1.4  IE-03 
8.56E-03 
8.79E-03 
-9.19E-04 
7.21E-04 
8.23E-04 
4.12E-03 
-4.93E-03 
-1.19E-03 
2.38E-03 
3.59E-03 
-2.89E-04 
-1.94E-04 
1.63E-06 
2.14E-06 
8.79E-04 
-4.13E-04 
2.42E-05 
5.03E-04 
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G  (n.m)  X  =  0°  X  -  -26.7° 

G  (10,10)  0  -1.39E-07 


